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Abstract 

The development and certification of airborne forward-looking wind shear detec- 
tion systems has required a hazard definition stated in terms of sensor observable 
wind field characteristics. This paper outlines the definition of the F-factor wind 
shear hazard index and an average F-factor quantity, F, calculated over a specified 
averaging interval, which may be used to judge an aircraft s potential performance 
loss due to a given wind shear field. A technique for estimating airplane energy 
changes during a wind shear encounter is presented and used to determine the 
wind shear intensity, as a function of the averaging interval, that presents a signifi- 
cant hazard to transport category airplanes. The wind shear hazard levels are com- 
pared to F values at various averaging intervals for four actual wind shear 
encounters. Results indicate that averaging intervals of about 1 kilometer could be 
used in a simple method to discern hazardous shears. 


1 Nomenclature 


D aircraft drag 

E aircraft total energy: kinetic energy + potential energy 

E' aircraft specific total energy 

e airspeed unit vector 

- a 

F F-factor 

F averaged F-factor 

g acceleration due to gravity 

h height above ground 



m 


mass of aircraft 


V ' airspeed 

W time rate of change of inertial wind vector W ] 

W 1 vertical component of the inertial wind (positive down) 

w aircraft weight = mg 

2 Introduction 

The development of airborne forward-looking wind shear detection systems and their certification 
standards has required that the wind shear hazard be defined. Characteristics of the required defi- 
nition include compatibility with existing reactive detection systems and crew training, that air- 
craft performance limitations be considered, and that the definition parameters be quantifiable by 
remote sensors. Existing reactive detection systems utilize the F-factor hazard index, suitably fil- 
tered for gust rejection (Reference 1). This paper outlines the definition of the F-factor wind 
shear hazard index, the derivation of the components contributing to the magnitude of the F-fac- 
tor, and the development of a quantity F which may be used to judge the aircraft performance loss 
potential of a given wind shear field. The F quantity may then be used to determine appropriate 
wind shear alerting thresholds following the definition of acceptable aircraft airspeed and height 
loss and a number of other important variables. 

3 Hazard Index Development 

3 . 1 The F-Factor Definition 

A full description of the F-factor, its implementation and its use, can be found in Reference 1. An 
abbreviated description is included herein. 

The total specific energy of an aircraft (relative to the airmass) is given by 

F y2 

- = E = ~+h (l) 

m 2 g 

and its rate of change by 
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Incorporating wind terms into the aircraft wind axis equations of motion, the rate of change of 
specific energy can also be shown to be 





T-D 
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where the F-factor wind shear index, F, is defined as 


(3) 


F = 
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— I -a 
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(4) 


The first term in the F-factor is the dot product of the derivative of the inertial wind vector and a 
unit vector along the aircraft airspeed vector divided by the acceleration of gravity. The second 
term is the local inertial vertical wind speed divided by aircraft airspeed. 


The F-factor describes the rate at which the wind field changes the energy of the aircraft. A nega- 
tive value of F is produced by a performance increasing wind shear such as an increasing head- 
wind or updraft (or a combination of the two), and a positive value of F by a performance 
decreasing shear such as a decreasing headwind or downdraft (or a combination of the two). An 
aircraft experiencing a wind shear with an F-factor of 0. 1 must produce a ( T—D)/w (or specific 
excess thrust) ratio of 0.1 in order to maintain level flight at constant speed. Conversely, a perfor- 
mance increasing shear of -0.1 F-factor requires a specific excess thrust reduction to maintain 
steady flight. Importantly for large magnitude wind shears, an F-factor value which exceeds the 
maximum specific excess thrust performance of an aircraft will force that aircraft to lose some 
combination of airspeed or altitude regardless of pilot input. For reference, a typical twin-engine 
commercial transport may have a maximum specific excess thrust performance of approximately 
0. 17, a typical three-engine aircraft 0.13, and a typical four-engine aircraft 0. 1 1 . All values are for 
maximum available thrust in a clean configuration at the aircraft's maximum gross weight. The 
values may vary depending on the aircraft type and configuration, however these numbers pro- 
vide a useful illustration of the calculation technique as well as a quantification of the relative 
effects of wind shear among two-, three-, and four-engine transport aircraft. 
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3.2 The F Index 


Reactive systems, which cannot predict the scale length of a wind shear, must rely on in situ mea- 
surements of the F-factor and apply gust rejection filters to minimize nuisance alerts. Forward- 
looking systems provide the capability to determine the spatial extent of a wind shear prior to the 
encounter. This additional information may be used to compute a hazard index that scales with 
the aircraft’s performance and energy margins. However the F-factor is proportional to the time 
rate of change of specific energy as shown in Equation (3). Therefore the aircraft’s temporal net 
energy loss or gain must be related to the measured spatial variations of F. 

Combining Equations (2) and (3) gives 
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T-D 
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-F = 


V V n 

a a 
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+ h 


(5) 


The measured quantities V , V a , and h, are functions of time and may be transformed from tem- 
poral variations to spatial variations calculated along the flight path, s. For level flight 


0 E _ dEds _ dE 
dt ds dt g ds 


from Equation (3) 


dE 

ds 


V a (T — D 
V V w 



8 

Airspeed and inertial groundspeed are related by (for level flight) 


( 6 ) 
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v s = 


( 8 ) 


By relating the relative magnitudes, it can be shown that for small flight path angles (y ) 


,-of » 


(9) 


8 ' 8 

At points within a wind shear the second term in this approximation may be significant, but on 
integrating over the length of a headwind/tailwind shear the approximation will hold. Therefore 
the approximation can be made 
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Integrating from a point in space, s Q , over some interval L gives 


A E'(s 0 ,L) = J 


s +L s„ + L 

O 
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ds - Fds 


( 11 ) 


From Equation (1) the total net energy change experienced by an aircraft passing through this 
region of shear (between s q and s q + L) is given by 


AE '= E s + L~ E s, 


( y 2 ^ 

— + h 
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( v 2 ^ 

, — + h . 
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or 


A (V 2 ) 

A E'(s o ,L) = — -^-+A h 

2 g 


where A (V a ) 


(V 2 ) - (V 2 ) . Combining Equations (11) and (13) gives 

O s + L a s 
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and 
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Defining the averaged integral quantity 
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gives 
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( 17 ) 


F(s 0 , L) represents an equivalent average F-factor over some interval L, starting at position s 0 . 


4 Applications 

4.1 Takeoff and Landing Performance 

The first application of Equation (17) will be to determine the typical wind shear penetration per- 
formance capabilities of two-, three-, and four-engine transport aircraft. Given specific energy 
input conditions, and an allowable energy loss, it is possible to determine the magnitude and the 
spatial extent of F along the flight path which an aircraft can withstand. 

The expression for F(s 0 , L) in Equation (17) may be used to determine the average wind shear 
intensities which, given a certain thrust response, will produce a specified performance loss 
expressed in terms of airspeed and altitude losses. In order to do this, assumptions must be made 
in the variation of {T- D) /w corresponding to a pilot recognizing the wind shear encounter, and 
responding to it by applying thrust. For takeoff it is assumed that the aircraft engines are develop- 
ing full power and the thrust is constant throughout the wind shear encounter. For an approach 
along a glideslope, the engines may be throttled back and the pilot’s response time to the onset of 
wind shear and the engines’ spool up time must be taken into account. A representative profile of 
{T- D) /w variation for wind shear encounters on landing is shown in Figure 1. 

Table 1 gives the parameters used in the calculation of F. The calculation begins at the beginning 
of the wind shear event, therefore s 0 = 0 (and in Figure 1, t = 0 represents the beginning of the 
event). The allowable airspeed loss in each case represents typical margins to stick shaker activa- 
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tion, and the allowable height loss is a postulated maximum acceptable worst case value. For a 
given value of the interval L the first term on the right hand side of Equation (17) is evaluated by 
determining the time elapsed since the beginning of the event using the initial airspeed 
(t = L/ V ) J , and integrating under the curve in Figure 1 from t = 0 up to t = L/ V . The other 
two terms in Equation (17) are calculated simply from the appropriate values given in Table 1. 

The results are shown in Figure 2. At intervals of less than about 500 m, extremely high F values 
are required for the limiting energy loss. At such short values of L, the curves for all aircraft types 
merge, suggesting that the individual aircraft performance characteristics do not greatly influence 
the energy losses. This result is strongest in the landing scenario, since the (T- D ) /w ratio is 
identical for all aircraft on the glideslope, and a finite time and distance is required to increase 
thrust. At interval values of 1 km, individual aircraft characteristics are apparent and the worst- 
case aircraft situation requires an F of 0.12 for limiting performance. Beyond 1 km, the curves 
rapidly asymptote to values close to the 1 km value. An exception is the two-engine landing sce- 
nario. At short intervals the engine spool-up is not completed while in the shear, lowering perfor- 
mance. At longer values the spool-up is completed well before exiting the shear, and the full 
performance capabilities of the aircraft are utilized. The engine spool-up factor is much less pro- 
nounced in the three- and four-engine cases, as the maximum performance of the airplane is 
essentially the same as the 1 km limiting F value. Note that the curves in Figure 2 represent 
results for the worst case low-altitude encounter. Much higher F values could be sustained given 
higher initial energy conditions. 

It is possible to relate the interval L to typical atmospheric phenomena in a general sense. High 
levels of F over short distances (up to 400 m) may be experienced as turbulence. Longer shear 
lengths (400 - 4000m) may be the results of microbursts, gust fronts, or other associated phenom- 
ena. Larger shear length (over 4000m) events may exist in the atmosphere but rarely at the levels 
shown to be hazardous to an aircraft. 

Note that for a given thrust profile, F values calculated are a function of initial and final airspeed 
and height conditions only. The energy loss an aircraft experiences due to a shear over a given 
distance is a function of the integral of the instantaneous shear values. No unique spatial distribu- 
tion of that shear is required in these calculations. Due to the non-linearity of Equation (5) (and 

1. In fact in a wind shear encounter the airspeed and groundspeed are changing and not necessarily equal. The fact that the initial 
airspeed is used in the space/time conversion illustrates that this is an approximation to a non-linear problem. 
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the consequent non-linearity of Equation (17)) several such simplifying assumptions had to be 
made. These assumptions, effectively linearizing the problem, allow a first approximation to the 
wind shear hazard threat definition to be made. A more complex analysis procedure than pre- 
sented above would be required to evaluate the full non-linear problem. Similar energy analysis 
techniques have been used in the past, however, and compared with non-piloted simulation results 
obtained from a point-mass, Boeing 737 model that incorporated realistic lift and drag changes 
with angle of attack (Reference 2). In that study the simulation results were consistent with the 
energy analysis. Results from the point mass Boeing 737 model have also compared well with 
full six degree of freedom, non-linear, piloted simulation study results (Reference 3). 

4.2 Wind Shear Case Studies 

The curves generated in Section 4.1 define an aircraft’s performance capability in the presence of 
energy losses and can be related directly to values of F from atmospheric measurements. The 
curves do not, however, provide guidance on the scale lengths of real wind shear events that 
should be used in threat determination. This section presents F as a function of averaging interval 
for four actual wind shear encounters by aircraft to provide insight on those scale lengths that can 
discriminate true threats. 

Calculation Method 

A discrete spatial series of F, {F n } , may be obtained from on-board flight data recordings (con- 
verted from time series) consisting of N points equally spaced by A r, (n = 1, 2, ..., N) . For a 
particular value of L (the averaging interval), discrete values of F(s o , L) are calculated using a 
discretized form of Equation (16) 

m + N 

ave 

{ Fj(L )} = A- X F J m= < N ~ N aJ < 18 > 

v ave 

j — m 

where N ave = L/ (A r) . The resulting series, {F. (L) } , is a moving average over N ave points, 
starting at the beginning of the series, F , , and ending at F N _ N . The value of interest (worst 

ave _ 

case shear) for a given averaging interval L is the maximum value of F in this series. This value 
may be compared to the aircraft performance curves of Figure 2. This method is now applied to 
four actual wind shear encounters. 
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Case.l 


A Boeing 737 aircraft on final approach to Denver Stapleton Airport which penetrated a 
microburst and descended below 100 feet AGL well short of the runway threshold prior to estab- 
lishing a positive rate of climb and recovery (Reference 4). The values of F calculated from the 
flight data recorder for this event is compared to a typical 2 engine landing performance curve in 
Figure 3. 

Cas£_2 

A Boeing 767 aircraft on final approach to Atlanta Hartsfield Airport which penetrated a 
microburst and descended to approximately 70 feet AGL short of the runway threshold prior to 
establishing a positive rate of climb and successful missed-approach. As in Case 1, the F curve 
calculated from the flight recorder data is compared with a typical 2 engine aircraft performance 
curve and is shown in Figure 4. 

Case 3 

A Lockheed L-1011 aircraft which crashed while on the approach to Dallas-Ft. Worth Airport in 
May 1985 (Reference 5). The aircraft entered a microburst while at approximately 750 feet AGL 
and was not able to recover. The derived curves are shown in Figure 5. 

Case A 

NASA’s Boeing 737 research aircraft which deliberately penetrated a microburst at approximately 
900 feet AGL near Orlando FL. in June 1991 (Reference 6). Figure 6 shows the F curves for this 
event as well as the appropriate 2 engine performance curve. A description of the NASA research 
flight program is found in Reference 7. 

Also included in the figures is an alert threshold boundary specified for reactive detection systems 
(Reference 8). As can be seen from the figures, all the above cases were of large enough magni- 
tude to warrant an alert. Case 3 was clearly more severe than the other three cases. Assuming a 
specific excess thrust performance capability for the L-1011 of 0.15, this microburst shear 
exceeded the aircraft’s ability to maintain level flight for nearly three kilometers, or nearly 40 sec- 
onds of flight time. 

As can be seen from these case studies, the F value of real events varies significantly with the 
averaging interval chosen. Every event shown had a significant F value at very short averaging 
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intervals, although not always in excess of the aircraft performance limit curve. Averaging inter- 
vals less than about 500 m would be a poor choice for discriminating threats due to the presence 
of turbulence and the poor comparison with aircraft performance. Reactive detection systems are 
designed to reject such small scale motion even with F values in excess of 0.2. Beyond averaging 
intervals of 1.5 to 3.5 km, F values fall below the reactive detection thresholds. Even large scale, 
severe events such as the Dallas - Fort Worth microburst do not sustain shears over greater aver- 
aging intervals than about 3 to 4 km. Averaging intervals greater than about 1.5 km to 2 km 
would also be a poor choice for a detection system, as threatening events of small spatial extent 
may produce very low F values at these averaging intervals. In each case study, the F value 
exceeded both the reactive threshold and the aircraft performance limit curve at averaging inter- 
vals near 1 km. 


5 Concluding Remarks 

The analysis presented above provides a suitable parameter, F, for quantification of the wind shear 
hazard by forward-looking sensors, as well as a method of establishing a hazard definition for sys- 
tems requirements and certification testing. Given aircraft performance characteristics, 
microburst encounter geometry, and a minimum allowable aircraft energy at microburst exit, a 
limiting F value may be found which will produce the specified exit energy. The energy analysis 
employed does not model all factors that are involved in a wind shear recovery, such as the pilot- 
ing technique used to trade altitude for airspeed or the variation of F-factor with varying paths, but 
does provide a first-order approximation to the survivability of a particular wind shear encounter. 
The analysis has both conservative and non-conservative factors. The results are conservative in 
that the performance-increasing shear, usually encountered prior to hazardous wind shears, is not 
modeled. This initial performance increase would provide additional energy for the recovery, and 
is particularly important in very small scale events, such as turbulence, where the aircraft crosses 
performance-decreasing regions and reenters performance-increasing regions before any signifi- 
cant flight path change can occur. This effect is not modeled and results in artificially low values 
of limiting F at short averaging intervals. The allowable energy loss values used in this analysis 
also effectively place the microburst in the worst-case location for landing and departure. Any 
other location would increase survivability. The results can be made even more conservative by 
assuming lower aircraft performance, longer engine spool-up time or pilot delay, or a higher 
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required minimum energy. The results are non-conservative in that adequate airplane perfor- 
mance for recovery does not guarantee survival, given variances in pilot technique. This effect 
may be minimized by specifying ample exit energy, particularly in the airspeed margins. 

Comparison of the performance limit curves with estimated F values from a variety of real events 
suggests that the best simple discriminator of a hazard is F taken at an averaging interval near 1 
kilometer. Lower averaging intervals will be dominated by turbulence, which is not a perfor- 
mance threat, and higher averaging intervals may result in very low F values for hazardous events 
that would activate an installed reactive detection system. In each case studied, the hazardous 
wind shear event exceeded the appropriate aircraft performance limit curve at scale lengths near 1 
kilometer. 

The analysis presented here may be used to derive forward-look sensor certification requirements 
and establish reasonable hazard thresholds, when combined with appropriate takeoff and landing 
configuration airplane performance characteristics, and minimum energy requirements acceptable 
to the FAA and the aviation community. 
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